基于BAP的岩体结构面粗糙度与峰值抗剪强度研究

【摘要】岩体是自然界中具有一定结构特征的地质体,其内部发育有各种各样具有一定方向、规模和形态的结构面,如层面、节理、断层、裂隙、不整合接触面等等。结构面的存在使得岩体不同于其他材料,物理力学性质存在不连续性、各向异性、非均一性等特点,同时也降低了岩体的完整性和整体强度。岩体往往不是因为本身材料破坏发生失稳,而是沿着岩体内部的软弱结构面发生破坏。甚至可以说,岩体的稳定性在一定程度上是由这些结构面(特别是控制性结构面)所决定的。因此,开展结构面抗剪强度力学参数的研究是评价岩体稳定性的关键所在。而作为影响岩体结构面抗剪强度的重要因素之一,粗糙度与抗剪强度之间具有明显的相关性,结构面越粗糙,抗剪强度越大。基于粗糙度的结构面峰值抗剪强度估算模型研究一直是岩体力学与工程地质领域的重要研究课题,估算模型的成败在于能否寻找出合理的指标来表征岩体结构面的真实粗糙程度。然而,目前关于结构面粗糙度的研究仍然存在以下几个问题：(1)采集装备落后,试验周期长,不能满足岩体稳定性的快速评价要求。对于发展迅速的岩体破坏过程,如高速远程滑坡,还有一些重大地质灾害抢险工程,时间就是生命,需要研究学者能够在短时间内给出评价结果,并预测受灾范围,提出合理的工程治理方案。而恰恰由于测量仪器本身问题,却体现不出估算模型的优势。(2)评价方法对采样方向、采样尺寸以及采样间距具有很大的依赖性,缺乏系统地开展粗糙度影响因素的研究及其相应合理的解决对策。(3)粗糙度评价指标仅仅依据结构面的几何信息,忽视了结构面粗糙度的力学特征,结构面粗糙度与结构面的受力状态、受力方向是密切相关的。本论文研究目的就是借助先进的采样手段,系统分析采样参数(采样方向、采样尺寸以及采样间距)对结构面粗糙度评价结果的影响,进而提出一种快速精确的粗糙度评价方法,并力求在此基础上建立结构面峰值抗剪强度估算模型。以重庆武隆鸡尾山高速远程滑坡研究区域的灰岩结构面为例,本论文主要开展以下几方面的研究工作：(1)天然岩体结构面粗糙度数据的快速采集。LS05-1岩体结构面为散落在堆积区内的灰岩层面,非原位产状,采用ILRIS-36D三维激光扫描仪,获取野外大尺寸结构面粗糙度点云数据,这部分数据主要用于岩体结构面粗糙度影响因素、基于光亮面积百分比(BAP)的粗糙度表征方法和结构面峰值抗剪强度估算模型研究工作中,原始点云数据扫描尺寸(ROI)为4.7×2.5m,采样间距为0.8mm,共采集扫描点约1800多万个。JWO2-1BL、JW02-2BL、JW03-1BL、JW03-2BL和JW03-3BR岩体结构面为采自于滑动面上的灰岩节理,运用VZ-400三维激光扫描仪,获取室内小尺寸结构面粗糙度点云数据,这部分数据主要服务于结构面峰值抗剪强度估算模型的研究工作,为数值模拟试验提供参数并验证计算结果,原始点云数据平均扫描尺寸(ROI)为0.181×0.362m,采样间距为1.0mm,每个结构面共采集65500左右个扫描点。对采集的原始点云数据进行预处理,主要包括点云数据解析、去噪和精简、几何模型建立三大部分,综合协调计算机运算能力、精度与处理时间等因素,最后LS05-1结构面精简成尺寸为4.70×1.14m,间距为0.01m,约54000个扫描点的点云数据；而JW02-1BL、JW02-2BL、JWO3-1BL、JW03-2BL和JW03-3BR结构面精简成平均尺寸为0.06×0.06m,间距保持不变仍为1.0mm,每块约3600个扫描点的点云数据,利用编写的Matlab程序,建立岩体结构面三维几何模型。(2)岩体结构面粗糙度影响因素研究。以LS05-1岩体结构面预处理后的点云数据为主,Weierstrass-Mandelbrot分形随机函数产生为辅,分别从二维与三维角度,分析采样方向、采样尺寸以及采样间距对结构面粗糙度评价结果影响。分析过程中做到单一变量原则,其余两因素保持固定。以15。为间隔,选取从0。到345。共24个采样方向,采用分形理论计算各个方向下的结构面二维剖面分形维数D,对于三维结构面,采用改进的赤平投影极点图法求出相应方向下的法向量数量NV。选取9种不同尺寸(长度或面积)的岩体结构面,基于分形理论求出相应尺度下的分形维数D。以0.01m为间隔,选取从0.01m到1.00m共100种采样间距,采用分形理论计算相应间距下的分形维数D。并且针对每一种影响因素,提出切实可行的解决方法。(3)基于光亮面积百分比(BAP)粗糙度表征方法的提出。从无充填低法向应力下的岩体结构面剪切破坏机制入手,以结构面中面向剪切方向的微小平面为研究对象,提出表征结构面粗糙度的新指标-BAP。基于Matlab软件平台,编写了光源模拟程序,并基于图像分割技术,给出BAP的定义和计算方法。对光源模拟与图像分割过程中涉及的两重要参数-光源入射角β和灰度阈值TL取值问题进行了探讨。通过描述结构面粗糙度的各向异性、尺寸效应以及间距效应,检验BAP法的适用性。(4)基于BAP的结构面峰值抗剪强度估算模型的建立。基于Matlab软件平台编写XYZ23DEC程序,简便快捷地建立起结构面三维复杂离散元模型。通过与室内直剪试验结果对比,来调整和验证数值试验所需的岩块、节理力学参数。在此参数基础上,以抗剪强度的各向异性、尺寸效应以及对法向应力敏感性等特征为参照点,验证数值模型的可靠性。设置合理的采样尺寸与采样间距常量,考虑结构面的粗糙度BAP、各向异性(剪切方向α)、岩块强度(JCS)、结构面基本摩擦角(φb)以及法向应力(σ)等参数,开展60例结构面直剪数值试验,统计分析模拟结果,建立粗糙度评价指标BAP与结构面峰值抗剪强度之间的关系,提出后者估算模型。(5)结构面峰值抗剪强度估算模型的应用。将该模型应用于重庆武隆鸡尾山高速远程滑坡岩体稳定性评价当中,结合研究区域具体的粗糙度BAP、法向应力σ、基本摩擦角φh以及壁面强度JCS参数,估算出滑动面抗剪力学参数(内摩擦角φ和黏聚力C)的变化区间。选择最小值、最大值、平均值三种组合,基于大型三维离散元软件3DEC,再现滑体初期破坏演化过程,分析其变形特征。在上述研究内容与方法的基础上,得出如下几点结论：(1)运用三维激光扫描技术可以实现岩体结构面粗糙度数据的快速采集目的。扫描时间与扫描间距、扫描范围直接相关。扫描间距越小,扫描范围越大,所花费的扫描时间会随之增加。除此之外,扫描时间还受到非技术因素的影响,如操作不当、电量供给、操作环境等等。扫描得到的原始点云数据需要对其进行预处理工作,点云数据的预处理是一个重复的过程,需要不断根据反馈进行修正与改善,预处理工作花费的时间很大程度上取决于工作人员的专业素质,操作越熟练所用时间越短。(2)采样方向对岩体结构面粗糙度评价影响分析显示：不同采样方向上的粗糙度之间存在较大差异,解决采样方向引入误差的最直接最有效方法,就是在数据采集过程中,确保采样方向与研究方向保持一致,研究方向包括滑动方向、剪切方向、渗流方向等等,这些研究方向需要在地质调查阶段进行判别,如滑动方向可根据滑动面上的擦痕和阶步等现象进行确定。采样尺寸对岩体结构面粗糙度评价影响分析显示：不同采样尺寸的结构面其粗糙度也各不相同,粗糙度会随着采样尺寸的增加而减小,当尺寸增加到有效尺寸(3.5m和3.5×0.84m)时,结构面粗糙度基本趋于平稳不再下降,通过增加采样尺寸的方法可以有效消除粗糙度的尺寸效应,但是受目前测量仪器的量程限制,也可以通过增加测量数量来弥补测量尺寸的不足,同样可以达到有效消除尺寸效应的目的,该思路在当前阶段具有较强的实用价值。采样间距对岩体结构面粗糙度评价影响分析显示：不同采样间距的结构面具有不同的粗糙度,整体上随着采样间距的减小,更多细节信息被获取,相应的粗糙度也越来越大,当采样间距小到有效间距时,粗糙度上升趋势变缓,趋于稳定状态,同时当采样间距过大时,大量信息被忽略,使得岩体结构面不再具备分形特征。有效采样问距与采样尺寸呈正比关系,建议采样间距与采样尺寸的比值不应高于0.05,在此条件下,可以有效降低采样间距对粗糙度评价结果的影响。(3)在无充填低法向应力作用下的岩体结构面剪切破坏过程中,面向剪切方向的微小平面主要发生剪切破坏,抗剪作用较大,结构面粗糙度表征方法的研究对象应集中在该接触部分。灰度阈值TL的选取至关重要,应结合结构面剪切破坏力学机制,遵循具体情况具体分析原则进行选取。模拟光源入射角β对BAP计算结果存在显著影响,BAP会随着入射角β的增大而增大,建议最佳入射角β取值范围为35。与70。之间。具体案例应用中,入射角β与灰度阈值TL的选取应该相互协调配合,给出合理组合。提出的BAP表征新方法可以很好地反映出岩体结构面粗糙度的各向异性、尺寸效应以及间距效应特性,该法是综合考虑结构面几何与力学因素下提出的粗糙度表征新方法,具有独特优势。(4)结合室内岩体结构面直剪试验展开的数值模拟结果显示：数值试验可以很好地体现出结构面抗剪强度的各向异性、尺寸效应以及与法向应力密切相关等特征。用于图像分割中的灰度阈值TL,其数值选取受法向应力σ的影响。以重庆武隆鸡尾山高速远程滑坡堆积区的LS05-1结构面为例,在光源入射角β=35°的前提下,法向应力所对应的合理灰度阈值TL如下所示,法向应力σ=0.5MPa时,灰度阈值TL=160;法向应力σ=1.0MPa时,灰度阈值TL=155;法向应力σ=1.5MPa时,灰度阈值TL=150;法向应力σ=2.0MPa时,灰度阈值TL=145;法向应力σ=2.5MPa时,灰度阈值TL=140。在大量数值模拟数据的基础上,经过多元统计分析,拟合出结构面峰值抗剪强度估算新模型,其形式如下,(5)通过27个案例估算出武隆鸡尾山高速远程滑坡滑动面9组黏聚力C和内摩擦角φ建议值,对比已出版文献数据,与各学者采用参数具有相似的统计规律,黏聚力C和内摩擦角φ的最小值组合为:φ=3.977°,C=0.443MPa,最大值组合为:φ=38.080°,C=1.917MPa,平均值组合为：(p=17.162°,C=1.094MPa,工况一(最小值组合)中提供的滑动面黏聚力C和内摩擦角φ比较符合实际情况。以工况一为基础,模拟滑坡破坏初期滑体的位移变化规律,数值模拟结果显示：滑体不同部位处的岩体具有不同的变形规律。滑体前缘和后缘块体具有较大的位移,前缘滑体沿T2裂缝走向滑动,后缘块体沿滑动面和T1裂缝交线方向滑动,滑体前缘块体由于临空面的存在较先发生破坏；相同位置处的滑体上部块体位移量比下部要大,上部块体先于下部块体发生破坏；滑体外侧(东侧)块体位移量要大于内侧位移量,尤其滑体后缘部分差异更为明显,不同步移动造成滑体破坏过程中发生旋转运动。通过与现场地质调查结果对比,证明了数值模拟结果的合理性,同时也证明了BAP结构面峰值抗剪强度估算模型得到了较为成功的应用。本论文主要创新成果有：(1)提出表征岩体结构面粗糙度的新指标BAP。它是一种三维评价方法,可以考虑粗糙度的各向异性、尺寸效应与间距效应等特征,并且结合了岩体结构面的剪切破坏机制；该法基于三维激光扫描试验获取的点云数据,所有计算过程实现了程序化,评价精度和速度得到大幅度提高。(2)基于三维粗糙度评价新指标BAP,建立岩体结构面峰值抗剪强度估算新模型。该模型是在合理的采样尺寸与采样间距下,考虑了与结构面相关法向应力σ、壁面强度JCS、基本摩擦角φb、粗糙度BAP以及剪切方向等因素而提出的；结构相对来说比较简单,具有很强的力学含义,模型中涉及到的统计参数是在大量数值试验结果基础上拟合而得,具有一定的可靠性；由该模型的实际工程应用显示,适用性较强,能够估算出合理的结构面抗剪力学参数。
【作者】葛云峰；
【导师】唐辉明；
【作者基本信息】中国地质大学，地质工程，2014，博士
【关键词】岩体结构面；粗糙度；峰值抗剪强度；三维激光扫描；数值分析；

【参考文献】
[1]章磊.基于RFCMOS工艺的2.4GHz正交频率合成器的研究与设计[D].西安电子科技大学,微电子学与固体电子学,2011,硕士.
[2]林晓静.山核桃质量控制中的统计方法应用研究[D].北京林业大学,统计学,2013,硕士.
[3]陈藻.农民工“半城市化”问题研究[D].西南财经大学,人口学,2013,博士.
[4]王丽.新型共聚聚乳酸/纳米β/磷酸三钙可吸收双层椎间融合器的研制及其降解性能研究[D].复旦大学,材料物理与化学,2012,硕士.
[5]刘化波.一种基于JPEG和HVS的水印方法[D].大连海事大学,应用数学,2004,硕士.
[6]桂凯旋.ZrB_2基纳米超高温陶瓷复合材料的制备及性能研究[D].哈尔滨工业大学,材料学,2013,硕士.
[7]唐敏.山区城市隧道节能减排技术研究[D].重庆交通大学,市政工程,2012,硕士.
[8]任中敏.不同水平的象棋棋手共同情境觉察的差异[D].河北师范大学,基础心理学,2004,硕士.
[9]肖波.液相脉冲放电制备金属陶瓷涂层的计算机仿真及粒子模拟[D].广东工业大学,材料加工工程,2013,硕士.
[10]杨文杰,胡明昊,唐振民,杨静宇.一种有效的基于约束的关联规则发现算法[J].南京理工大学学报(自然科学版),2005,01:109-112.
[11]王丹阳.不同性别化女大学生的客体关系表征水平[D].河北师范大学,应用心理,2014,硕士.
[12]王红梅,王蔺,王晓川.点击反应在功能材料中的应用[J].材料导报,2014,15:130-135.
[13]胡倩.明代宗室的文化成就研究[D].湖南师范大学,历史文献学,2013,硕士.
[14]王峭.敏化剂改性纳米TiO_2的光催化性能研究[D].兰州理工大学,材料学,2014,硕士.
[15]肖登元.论中国山水画与现代设计相结合的逻辑[D].西安美术学院,美术学,2014,硕士.
[16]刘欢欢.大功率DC-DC变换器串并联技术的研究[D].西南交通大学,电力电子与电力传动,2013,硕士.
[17]DilanthaThushara（苏贝）.贵金属高效非均相催化剂Pd-BTP/SiO_2和Pd-石墨烯在Suzuki和Heck反应中的应用[D].浙江大学,2014.
[18]王汀,江金财,俞志勇,徐骁.2013中国器官移植大会会议纪要[J].中华移植杂志(电子版),2014,01:45-49.
[19]肖锋,吴燕冈,孟令顺.位场数据处理中小子域滤波的改进[J].石油地球物理勘探,2010,01:136-139+164+175.
[20]汪连海.甘草快速繁殖技术及影响因素的研究[D].甘肃农业大学,作物栽培学与耕作学,2013,硕士.
[21]刘顺喜.高光谱遥感土地利用信息提取技术研究[D].北京林业大学,2005.
[22]伍晓榕.零件工艺失效模式关联建模与工艺方案绿色性分析关键技术研究[D].浙江大学,机械设计及理论,2014,博士.
[23]付岳.白藜芦醇抑制辐射诱导的NLRP3炎症蛋白复合体激活的机理研究[D].北京协和医学院,放射医学,2014,博士.
[24]高云芳,宋云龙,任冬雷,王艳平.纳米Pb(PbO)/活性炭的制备及其在H_2SO_4电解液中的电化学行为[J].无机材料学报,2013,12:1301-1306.
[25]俞佳玉.江苏省大中型工业企业技术创新能力评价研究[D].苏州大学,企业管理,2013,硕士.
[26]姚元.低成本战略更能提升造纸企业的竞争力吗？[D].北京交通大学,2014.
[27]任尊松,徐宇工,王璐雷,邱英政.强侧风对高速列车运行安全性影响研究[J].铁道学报,2006,06:46-50.
[28]赵倩.挂于腰间的金色文化甘肃保安族腰刀多元一体设计文化探析[D].重庆大学,设计学,2014,硕士.
[29]刘媛媛.蛴螬抗肿瘤活性化学成分的筛选与分离[D].第四军医大学,中药学（专业学位）,2014,硕士.
[30]史方舟.三维GIS输电线路管理系统的研究与应用[D].东北师范大学,地图学与地理信息系统,2012,硕士.
[31]张礼才.网络突发事件政府信息发布研究[D].华中科技大学,公共安全预警与应急管理,2013,硕士.
[32]王浙君.应采用全程式监理[J].中国计算机用户.2005(31)
[33]吴德昌.铁路管控系统QoS通信解决方案[D].西南交通大学,电力系统及其自动化,2003,硕士.
[34]闫晶晶.师范生免费教育政策的报纸媒介舆论环境研究[D].首都师范大学,教育经济与管理,2013,硕士.
[35]罗英花.基于价值链分析的DH公司战略转型研究[D].北京交通大学,工商管理（专业学位）,2013,硕士.
[36]王锦标,朱因心.图形组态软件中的递归式回路识别方法[J].冶金自动化,2003,01:35-38.
[37]乔姗.基于知识审计的国际货代业务流程建模与分析[D].大连海事大学,工业工程,2014,硕士.
[38]李金国.基于应用的高速网络入侵检测系统的研究与实现[D].湖南师范大学,2004.
[39]刘晓亮,类延达,郭宝春,贾德民.离子液体改性炭黑增强丁腈橡胶的研究[J].弹性体,2009,05:1-5.
[40]何天明.WAC-Ⅱ工业微机在杭钢一号高炉无料钟炉顶自动控制中的应用[J].冶金自动化,1993,06:44-45.
[41]徐海霞.汉语国际教育专业本科课程设置中的外语课程研究[D].华中师范大学,汉语国际教育,2014,硕士.
[42]刘海波.邻近基坑开挖对千厮门大桥引桥安全性影响评估[D].重庆大学,土木工程,2014,硕士.
[43]冯博.太阳能电池片表面裂纹检测算法的研究与实现[D].吉林大学,计算机应用技术,2014,硕士.
[44]王亚红.小额贷款公司发展中的法律困境及其出路[J].上海市经济管理干部学院学报,2014,04:43-50.
[45]金春红.中国商业银行经营效率实证研究[D].辽宁大学,金融学,2014,博士.
[46]王双.基于羧酸配体在离子液体中构筑的配位聚合物的合成、结构与性质研究[D].吉林大学,无机化学,2014,硕士.
[47]邵鸿玲.基于数据仓库技术在地铁故障系统中的应用[D].江苏科技大学,计算机应用技术,2013,硕士.
[48]王荣荣.孔子与柏拉图财富伦理思想的比较研究[D].宁夏大学,外国哲学,2014,硕士.
[49]吕祥.UCM冷连轧机薄带钢轧制板形控制的研究及有限元仿真[D].燕山大学,机械设计及理论,2013,硕士.
[50]陈发有.防屈曲支撑在有局部收进高层钢结构中的减震性能分析[D].兰州理工大学,结构工程,2014,硕士.